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Relic neutralinos and the two dark matter candidate events 
of the CDMS II experiment 
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The CDMS Collaboration has presented its results for the final exposure of the CDMS II ex- 
periment and reports that two candidate events for dark matter would survive after application of 
the various discrimination and subtraction procedures inherent in their analysis. We show that a 
population of relic neutralinos, which was already proved to fit the DAMA/LIBRA data on the 
annual modulation effect, could naturally also explain the two candidate CDMS II events, if these 
are actually due to a dark matter signal. 

PACS numbers: 95.35.+d,11.30.Pb,12.60. Jv,95.30.Cq 



The search for a sign from Dark Matter (DM) involves 
direct detection, consisting in the measurement of the 
effects induced by the feeble interaction of the DM par- 
ticles with the material of a low-background set-up, and 
indirect measurements. These concern many possible 
signals, ranging from neutrinos to charged cosmic rays 
(positrons, antiprotons, antideuterons), to gamma rays, 
to a radio signal and even to effects induced on the cosmic 
microwave background. 

In Rcf. [l[ we have shown that the annual modulation 
effect at a 8.2 a C.L., obtained by the DAMA/Nal and 
DAMA/LIBRA experiments (with a total exposure of 
0.82 ton yr) is very well fitted by relic neutralinos 
in an effective Minimal Supcrsymmetric Standard Model 
(effMSSM) at the electroweak scale defined in terms of a 
limited number of parameters. We recall that the effect 
measured by the DAMA Collaboration is the first and 
up-to-now unique evidence for a signal compatible with a 
typical signature (annual modulation) expected for dark 
matter particles [3[. 

Other experiments of WIMP direct detection do not 
currently have the capability of measuring the annual 
modulation effect and usually provide upper bounds 
for the expected signals These limits are obtained 
through complex procedures for discriminating electro- 
magnetic signals from recoil events and through delicate 
subtractions meant to separate putative WIMP signals 
from neutron-induced events. A major critical point in 
these experiments and related analyses is that the very 
signature (the annual modulation) of the searched signal 
cannot be employed in extracting the authentic events. 
Other potential difficulties are related to stability fea- 
tures and determination of the threshold and of the en- 
ergy scale. 

In Ref. [l[ it was also pointed out that the inclusion 



of these upper bounds, taken at their face value, would 
anyway allow a compatibility with the DAMA data for a 
range in the WIMP (neutralino) mass around 7-10 GeV. 
A similar result has been obtained also in Ref. 

The CDMS Collaboration has now presented its results 
for the final exposure of the CDMS II experiment . In 
that paper it is reported that two candidate events for 
DM would survive after application of various discrimina- 
tion and subtraction procedures, though a probability of 
23% exists that they are of a more prosaic origin. These 
2 events have recoil energies of 12.3 keV and 15.5 keV. 

If one assumes that the two candidate events are due 
to a WIMP particle with a coherent interaction with nu- 
clei, taking into account the CDMS total exposure, one 
can derive that the relevant 90% C.L. region in the plane 
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(up to a WIMP mass of 100 GcV) is the 



one displayed in Fig. 1 (m^ denotes the WIMP mass, 
CTg™!^'™"' is the WIMP-nucleon coherent cross section 
and ^ is the WIMP local fractional density) . In this Fig- 
ure, due to the low statistics, we have adopted the simple 
criterion to require n = 2 WIMP events (0.6 < n < 4.7 
at 90% C.L. for a Poissonian distribution) in the to- 
tal range of the recoil energy En observed by CDMS, 
10 keV < Eb, < 100 kcV. This is sufficient to capture 
the main features of the allowed region. This is also true 
for < 8.5 — 9 GeV where, depending on the values of 
the escape velocity in the Galaxy and on the rotational 
velocity of the Solar System, the event with Er = 15.5 
keV could in principle not be ascribed to a WIMP. In 
fact in this case the region shown in Fig. 1 overlaps with 
the one (not shown in Figure [T]) obtained by requiring 
only one WIMP event (0.11 < n < 3.44 at 90% C.L.) 
for 10 kcV < Eji < 12.3 kcV: the region obtained with 
this criterion has the upper boundary about 25% smaller 
than the one shown in Fig. [T] and the lower boundary 
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FIG. 1: £(J_^"7';°"> as a function of the WIMP mass. The 
(green) shaded regions denote the DAMA/LIBRA [31 annual 
modulation regions, under the hypothesis that the effect is 
due to a WIMP with a coherent interaction with nuclei; the 
region delimitated by the solid line refers to the case where 
the channeling effect is not included, the one with a dashed 
contour to the case where the channeling effect is included 
The (violet) band displays the region related to the two 
CDMS candidate events, obtained from the total rate in the 
whole energy window. The scatter plot represents supersym- 
metric configurations calculated with the model summarized 
in the Appendix, at a fixed representative set of values for 
the hadronic quantities. The (red) crosses denote configura- 
tions with a neutralino relic abundance which matches the 
WMAP cold dark matter amount (0.098 < O.^h'^ < 0.122), 
while the (blue) dots refer to configurations where the neu- 
tralino is subdominant {Q^h^ < 0.098). The region covered 
by a (blue) slant hatching denotes the extension of the scatter 
plot upwards and downwards, when the hadronic uncertain- 
ties in the scattering coherent cross-section are included. 

reduced by about a factor of 5. 

In Figure [1] also the annual modulation regions of the 
DAMA Collaboration are shown, with and without inclu- 
sion of the channeling effect 0. The exact modeling of 
channeling is still under study, then one expects that the 
actual physical situation is comprised within the two re- 
gions represented in the figure. As a reference model for 
the WIMP halo distribution, a cored-isothermal sphere 
is employed with the following parameters: local value of 
the rotational velocity vq ~ 220 km s^^, escape velocity 
Vesc = 650 km s~^ and total non-baryonic dark-matter 
density po = 0.34 GeV cm~^. Obviously, the DM halo 
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FIG. 2: The same as in Fig. [Jl except that the (yellow) shaded 
regions compatible with the CDMS II candidate events are 
obtained by a maximal likelihood method applied to the dif- 
ferential energy recoil rate, under the hypothesis of negligible 
background. The contours refer to (from the internal to the 
external one) 68%, 90% and 95% C.L. 

distribution could be quite different this would in- 
duce a shift of the actual position of the regions and 
bounds, as discussed and shown e.g. in Ref. jl|. 

In Fig. 1 we also display the scatter plot represent- 
ing the supersymmetric configurations calculated within 
a realization of the Minimal Supersymmetric Standard 
Model where gaugino unification is relaxed [l|, For 
convenience, the model is summarized in the Appendix. 
The scatter plot refers to a fixed representative set of val- 
ues for the hadronic quantities involved in the neutralino- 
nuclcon cross sections [l[ . The region covered by a (blue) 
slant hatching denotes the extension of the scatter plot 
upwards and downwards, when the hadronic uncertain- 
ties extensively discussed in Ref. [l| are included. 

From Fig. [T] and the previous discussion about the 
CDMS region one finds that the putative CDMS events 
are compatible simultaneously with the DAMA/LIBRA 
data and the theoretical evaluations in the mass range 
8-12 GeV. It is worthwile to point out that at such low 
masses the expected recoil spectrum depends on the high 
velocity tail of the velocity distribution, which is sensitive 
to the details of the astrophysical model and in particular 
to the escape velocity. Other possibilities for the mod- 
eling of the velocity distribution have been discussed in 
Q . We also stress that the explanation - in terms of the 
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FIG. 3: The same as in Fig. [2l but under the hypothesis of 
a backgound contribution as in Ref. [l^, normaUzed to 0.8 
events in the whole energy window of CDMS II. The contours 
refer to (from the internal to the external one) 68% and 85% 
C.L. 



relic neutralinos in the supersymmetric model discussed 
here - of the annual-modulation data alone extends over 
a much wider range; for instance, for the case of a WIMP 
halo distribution given by a cored-isothermal sphere with 
the parameters mentioned above, this extended range can 
be simply read from Fig. 1 to be 6 GeV < rriy^ < 60 GeV. 
These light neutralinos can also be complementarily in- 
vestigated by indirect means, such as cosmic antiprotons 
0, and antideuterons , signals at neutrino tele- 

scopes 
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and, most notably, can be searched for at the 
Large Hadron Collider [l^. Astrophysical bounds aris- 



ing from multi-wavelength analyses IJ], which may be 
strong depending on assumptions on the DM distribu- 
tion and on astrophysical properties, like those related 
to cosmic-rays propagation and energy losses, do not 
markedly constrain the supersymmetric configurations of 
Fig. [TJ especially when astrophysical uncertainties are 
properly taken into account. 

Our previous analysis was based only on the total rate 
taken over the whole recoil energy range observed by 
CDMS II, without using any spectral information. This 
is motivated by the very low statistics (2 events), which 
makes very critical (and to some extent not fully justi- 
fied) a statistical analysis of the energy spectrum. How- 
ever, forcing somewhat the situation, one can wonder 
what would produce an analysis in terms of the energy 



spectrum. In Fig. [2] and [3] we therefore show the re- 
gions compatible with the 2 CDMS II events at 12.3 keV 
and 15.5 keV, taking into account the spectral behaviour 
of the theoretical recoil rate. In the determination of 
the allowed regions we have adopted a maximal likel- 
hood analysis [l7[ . Fig. [5] shows the case of a negligi- 
ble background contribution, and the contours refer to 
a confidence level of 68%, 90% and 95%, from the in- 
nermost to the outermost. Fig. [3] instead refers to the 
presence of a back grou nd contribution, which we have 
modeled as in Ref. [18[, i.e. with an energy dependence 
dN/dE = -0.00295 + 0.463/S normalized to the total 
number of event of 0.8, to conform to an estimate of the 
background contribution of 0.8 ± O.l(stat) ± 0.2(syst)[^. 
In the case of Fig. [31 the contours refer to 68% and 85% 
C.L., and they evolve into an open region [i.e. into an 
upper bound) at the 90% C.L., a result compatible to 
the one obtained in Ref. [l8|, where a slightly different 
statistical analysis is adopted. Fig. [2] and [3| are quite 
compatible with the results of Fig. [1] obtained by using 
the total counting number, and reinforce our conclusions 
of compatibility between DAMA and CDMS II for light 
WIMPS, and between these experimental results and our 
SUSY models with light neutralino dark matter. 

In conclusion, in this note we have considered the two 
events which, in the analysis of CDMS II, seem to survive 
after the various discrimination and subtraction proce- 
dures. We have shown that, should these events be 
due to WIMP-nucleus coherent interactions, this result 
would be compatible both with the annual-modulation 
signal previously reported by the DAMA Collaboration 
and with an interpretation in terms of relic light neu- 
tralinos. This conclusion is not affected by other upper 
bounds of direct dark matter detection. In particular. 



the XENON upper bound [15| suffers from large uncer- 
taities due to conflicting determinations of the scintilla- 
tion efficiency at low nuclear recoils (as shown in Fig. 
12 of [l3|)- Calculations performed in Ref. [l^ (though 
with a threshold energy somewhat smaller than the one 
of XENONIO) indicate that at a WIMP mass of order 10 
GeV the bound of Ref. [l^l should be relaxed by more 
than an order of magnitude. One should furthermore 
note that in XENONIO the energy scale is particularly 
uncertain due to a calibration at an energy much higher 
than the declared threshold energy. 



I. APPENDIX: THE SUPERSYMMETRIC 
MODEL 

The supersymmetric scheme we employ in the present 
paper is the one described in Ref. [l[ as an effective 
Minimal Supersymmetric Standard Model (cffMSSM) at 
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the electroweak scale, with the following independent pa- 
rameters: Ml, M2, M3, /X, tan/3, m^, m,, and A. No- 
tations are as follows: Mi, M2 and M3 are the U(l), 
SU(2) and SU(3) gaugino masses (these parameters are 
taken here to be positive), /x is the Higgs mixing mass 
parameter, tan/3 the ratio of the two Higgs v.c.v.'s, 
the mass of the CP-odd neutral Higgs boson, rriq is a 
squark soft-mass common to all squarks, mj is a slepton 
soft-mass common to all sleptons, and A is a common 
dimensionlcss trilincar parameter for the third family, 
A-j^ = Aj: = Arriq and Af = Amj (the trilinear param- 
eters for the other families being set equal to zero). In 
this model no gaugino mass unification at a Grand Uni- 
fied scale is assumed, whence light neutralinos arise. 

The parameter space of this model is bounded by a 
large host of experimental data: invisible Z decay (for 
decay into light neutralinos), direct searches of super- 
symmetric particles and higgs bosons at LEP and Teva- 
tron, supersymmetric constributions to rare processes: 
BR{b — > s + 7), BR{B^ — > + /i+), measurements of 
the muon anomalous magnetic moment = (p^ — 2)/2. 
Other details about the model and the relevant con- 
straints can be found in Ref. 



II. NOTE ADDED 

After submission of the present paper, the CoGeNT 
Collaboration [l^ has presented the results of a search 



for light-mass DM particles, where an irreducible excess 
of bulk-like events below an energy of 3 keV is observed. 
As discussed in Ref. [l^ , should this population of events 
be due to WIMP interactions with the detector, these 
would entail a WIMP mass of 7-12 GeV with a cross 



section m. 



>. (nuclcon) 
^ ^scalar 



(3 - 10) X 10- 



thus in a 



region in agreement with predictions of our model, with 
the DAMA/LIBRA and CDMS II resuhs. 
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